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Abstract 

In the context of some deformed canonical commutation relations leading to 
isotropic nonzero minimal uncertainties in the position coordinates, a Dirac equation 
is exactly solved for the first time, namely that corresponding to the Dirac oscillator. 
Supersymmetric quantum mechanical and shape-invariance methods are used to de- 
rive both the energy spectrum and wavefunctions in the momentum representation. 
As for the conventional Dirac oscillator, there are neither negative-energy states for 
E = —1, nor symmetry between the I = j — ^ and I = j + ^ cases, both features being 
connected with supersymmetry or, equivalently, the lo —lo transformation. In con- 
trast with the conventional case, however, the energy spectrum does not present any 
degeneracy pattern apart from that associated with the rotational symmetry. More 
unexpectedly, deformation leads to a difference in behaviour between the / = j — ^ 
states corresponding to small, intermediate and very large j values in the sense that 
only for the first ones supersymmetry remains unbroken, while for the second ones no 
bound state does exist. 
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1 Introduction 



Many years ago, Ito et al [1] considered a Dirac equation in which the momentum P is 
replaced by P — \muj(5X with X being the position vector, m the mass of the particle and 
ui the oscillator frequency. It was later shown by Cook [2] to present unusual accidental 
degeneracies in its spectrum, which were discussed from a supersymmetric viewpoint by Ui 
and Takeda [3] and by Balantekin [4]. 

The interest in the problem was revived by Moshinsky and Szczcpaniak [5] , who gave it 
the name Dirac oscillator (DO) because in the nonrelativistic limit it becomes a harmonic 
oscillator with a very strong spin-orbit coupling term. Since then, the DO has aroused a 
lot of interest both because it provides one of the few examples of Dirac equation exact 
solvability and because of its numerous physical applications. 

As a relativistic quantum mechanical problem, the DO has been studied from many 
viewpoints, including covariance properties [6], complete energy spectrum and correspond- 
ing wavefunctions [7], symmetry Lie algebra [8], shift operators [9], hidden supcrsymme- 
try [7, 10, 11], conformal invariance properties [12], as well as completeness of wavefunc- 
tions [13]. 

Relativistic many-body problems with DO interactions have been extensively studied 
with special emphasis on the mass spectra of mesons (quark-antiquark systems) and baryons 
(three-quark systems) (see, e.g., [14, 15] and references quoted therein). The dynamics 
of wavepackets in a DO has been determined [16, 17] and a relation with the Jaynes- 
Cummings model established [17]. The 2-1-1 space time has also been shown [18] to be an 
interesting framework for discussing the DO in connection with new phenomena (such as the 
quantum Hall effect and fractional statistics) in condensed matter physics. Thermodynamic 
properties [19] of the DO in 1 + 1 space time [20] have been mentioned to be relevant to 
studies on quark-gluon plasma models. 

Various extensions of the DO have been considered in the literature. Let us mention 
some alternative approach based on a scalar coupling [21], as well as some generalizations 
to arbitrary spin [22] or to quasi-exactly solvable systems [23]. 

In the present paper, we plan to solve the DO problem in an entirely different context, 
namely that of theories assuming a modified Heisenberg uncertainty relation leading to 
a nonzero minimal uncertainty in position. Such theories have been stimulated by several 
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independent lines of investigation in string theory and quantum gravity suggesting the exis- 
tence of a finite lower bound to the possible resolution of length AXq (see, e.g., [24, 25, 26]). 
In particular, we shall be interested in a realization of this modified uncertainty relation 
through small quadratic corrections to the canonical commutation relations [27, 28]. An- 
other interesting aspect of these corrections is that they can also provide an effective de- 
scription of non-pointlike particles, such as quasiparticles and various collective excitations 
in solids, or composite particles, such as nucleons and nuclei [28]. 

Only a few nonrelativistic quantum mechanical problems have been considered in such 
a context so far. An exact solution to the one-dimensional harmonic oscillator problem 
has been provided by solving the corresponding deformed Schrodinger equation in momen- 
tum representation [29]. This approach has been extended to D dimensions [30]. Some 
perturbative [31] or partial [32] results have also been obtained for the hydrogen atom. 

The oscillator results, both in one and D dimensions, have been rederived [33, 34] by us- 
ing powerful combined techniques of supersymmetric quantum mechanics (SUSYQM) [35, 
36] and shape invariance (SI) under parameter translation [37, 38]. Such an approach 
is known [36, 39] to be a reformulation of the factorization method dating back to 
Schrodinger [40] and Infeld and Hull [41]. The combined formalism of SUSYQM and SI 
(this time under parameter scaling [42, 43]) has proved especially useful [33, 34] when one 
assumes nonzero minimal uncertainties in both position and momentum [44] , in which case 
neither position nor momentum representation can be used. 

Here we will avail ourselves of similar kinds of techniques to provide an exact solution to 
the DO problem with nonzero minimal uncertainty in position. In section 2, this problem 
is reviewed in momentum representation. In section 3, the equations for the large and 
small component radial momentum wavefunctions are derived and the corresponding SUSY 
partners obtained. The existence of physically acceptable wavefunctions corresponding to 
a vanishing energy for one of them is discussed in section 4. The DO energy spectrum and 
wavefunctions are then determined in sections 5 and 6, respectively. The final results are 
collected in section 7, while section 8 contains the conclusion. 

2 Dirac oscillator in momentum representation 

In a system of units wherein h — c — m — 1, the (stationary) DO equation is given by [5] 

Hi/j^Ei/j H ^a-{P- iuXp) + /3 (2.1) 
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where a hat denotes a 4 x 4 matrix, 

and (7j, i — 1, 2, 3, are the Pauh spin matrices. Here we assume that the position and 
momentum components X^, P^, i — 1 2, 3, satisfy modified commutation relations of the 
type [28, 29, 30] 

[Pi,Pj] = (2.3) 
[Xi,Xj] = -i[2/3-/3' + (2/3 + /3')/?P']e,,fcLfe 

where there is a summation over dummy indices, 

^ TTPP^''^'^'^' ^ = 1, 2, 3 (2.4) 

are the components of the orbital angular momentum, satisfying the usual commutation 
relations 

[Li, Xj] = ieijkXk [Li, Pj] = ieijkPk (2.5) 

and (3, (3' are two nonnegative, very small deforming parameters. In the DO problem, we 
shall assume that the two lengths and which are the square roots of the deforming 
parameters, are very small as compared with the oscillator characteristic length or, 
equivalently, that /3a; <S 1 and P'u <^ 1. 

The commutation relations (2.3) generalize to three dimensions (see, e.g., [28] for a 
derivation in d dimensions) the modified commutation relation between position and mo- 
mentum in one dimension 

[X,P]=i{l + (3P^). (2.6) 

The latter has been proved [29] to lead (in those states for which (P) = 0) to the uncertainty 
relation 

AX>^(^+/3AP) (2.7) 

appearing in perturbative string theory and in line with the proposed UV-IR mixing [24, 
25, 26]. Equation (2.7) implies a nonzero lower bound for AX: AX > AXq = y/p. 

For the commutation relations (2.3), the same mechanism as in one dimension yields 
from the corresponding uncertainty relation AXjAPj > ||([Xi,Pj])| an isotropic nonzero 



minimal uncertainty AXq = AXoi, i = 1, 2, 3, equal to AXq = \/3/5 + (3' when one restricts 
oneself to those states for which (Pj) = 0, i = 1, 2, 3, and AP, is independent of i [28]. 

Since, from (2.3), it follows that the momentum components remain simultaneously 
diagonalizable, we can work in the momentum representation, wherein Pj, and Lj are 
realized as 



Pi = Pi 
X, = i 



d d 
dpi dpj 



{21 



_ . d _ . d 

J^i — ^^ijkT{ — Pk — ~'^^ijkPjT\ — • 



Here 7 is an arbitrary constant, which does not appear in the commutation relations (2.3) 
and only affects the weight function in the scalar product in momentum space 



(C1C) 



[/(P)] 



where 



f{p) = l + (3,p^ (3o = f3 + i3' a = 



(2.9) 



(2.10) 



^1 



^2 



On separating the wavefunction tp = 
ponents, the DO equation can be written as two coupled equations 

where 

= CT ■ P ±iuj(T ■ X. 



in (2.1) into large ipi and small ip2 com- 



(2.11) 
(2.12) 

(2.13) 



Applying B+ (resp. B-) to (2.12) (resp. (2.11)) and using (2.11) (resp. (2.12), we get the 
following factorized equation for the large component ipi (resp. small component 1^2) 



B+B-ipi 



{E^ - l)V^i 



(2.14) 

S-S+V^2 = (^'-1)V^2. (2.15) 

As shown in appendix 1, the operators B^ and B~ can be written in the momentum 
representation (2.8) as 



B-" 



I P 

' p 



B = tUCTp < — h 



5 



P J dp 



P . 



d 



+ PP g^+lP 



(2.16) 
(2.17) 



where 



(7p = ^ (2.18) 



is such that 

al = I. (2.19) 

Some further simphfications are achieved by transforming the large and small compo- 
nents according to 

V^i = -r°/Vi ^ = l,2 (2.20) 
P 

where / and a are defined in (2.10). As a consequence, equations (2.11) and (2.12) become 

i3+02 = (^-1)01 (2.21) 
B-(l)i = (£^ + 1)02 (2.22) 

where denote the transformed operators i?^, which can be written as 

B+ = u;b+ap^u;apb+ (2.23) 
B~ = ouapb- = ub-ap. (2.24) 

Here the operators 6=*= are directly obtained from (2.10), (2.16), (2.17) and (2.20) as 

,± <9 , /I . , (t l + 1 



while the alternative forms 



r± . a , /I , . . cr- L + 1 



b^^:ff—+l- + (3(cT-L + l)\p + — (2.26) 

op \uj y p 

result from (2.25) and the relation 

{ap,(T-L + l}^0 (2.27) 

proved in appendix 1. Prom (2.19), (2.23) and (2.24), it follows that equations (2.21) and 
(2.22) can be rewritten as 

c^6+02 = (^-1)01 (2.28) 

u;b-(j)i = (E + l)4>2 (2.29) 

where 

02 = ap(j)2. (2.30) 



3 Large and small component radial momentum wave- 
functions and SUSY partners 



In section 2, we have demonstrated that the solution to the DO equation (2.1) with modified 
commutation relations (2.3) amounts to that of the system of coupled equations (2.28) and 
(2.29). Since the operators appearing on the left-hand side of such equations have the 
property of commuting with the total angular momentum J = L + S, where S = |cr, 
as well as with and S'^, we may look for solutions 0i and 02 that are simultaneous 
eigenfunctions of L'^, S^, and J3, corresponding to the eigenvalues l{l + 1), |, j(j + 1) 
and m, respectively. Instead of Z, we may use s, defined by I — j — s and taking the values 
±|. Then 0i and 02 can be expressed in momentum spherical coordinates p, 9, (p and spin 
variable ^ as 

01 = (l>l;s,j,m{P,0,ip,^) ^ Rl;s,j{p)ys,j,m{0,'P,O (3-1) 
h = 4>2;s,j,m(p, 9,ip,^)^ R2;s,j (p)ys,j,m(9, ifi, (3.2) 

where 

ys,j,m{^, V, = 51 (i -s n,^a\j m)Yj_,^f,{e, ip)xa{0 (3-3) 

fl,<T 

is a spin spherical harmonic [45] and Ri-sjip), R2;s,j{p) (or, in short, Riip), R2{p)) are 
radial wavef unctions. Note that in (3.3), Xo-(0> where cr = ±|, denotes a spinor. 

We note that 0i and 02 have the same angular-spin wavefunction. The same is not true 
for 01 and 02 since, from (2.19) and (2.30), the latter is given by 

02 = (^p02 = R2;s,j{p)crpys,j,m{0, (fi, = -R2;s,j{p)y-s,j,m{0, (fi, (3-4) 

where, in the last step, we used a well-known property of spin spherical harmonics [46]. 
Hence we may write 

02 = (j)2;-s,j,m{p, 0, (fi, = i?2;-sj (p)3^-s,i,m(^, (3-5) 

where 

R2--s,j{p) = -R2;s,j{p)- (3.6) 

In other words, 0i and 02 are characterized by the same / value {I = j — s), whereas 02 
corresponds to /' = j -|- s. 
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Inserting (3.1) and (3.2) in (2.28) and (2.29) and taking the property 

{(T-L + i)y,j,^ = ( j2 -L^-s^ + i)y,j,m = s{2j + i)y,j,^ (3.7) 

into account, we obtain a system of coupled radial equations 

u;b+R2 ^{E- l)i?i (3.8) 
ujb-Ri = {E+1)R2. (3.9) 

Here bp denote the radial parts of the operators 6=*= 

b^ = b^{g,k) = Tf^ + 9P-^ (3.10) 

where we have introduced the notations 

g^--f3s{2j + l) k^s{2j + l). (3.11) 

According to (2.20), (3.1) and (3.2), the scalar product (2.9) in momentum space gives 
rise to a scalar product 

{R'\R) = f ^^R'*ip)Rip) (3.12) 

in radial momentum space. It is easily checked that with respect to the latter the operators 
bp and b~ , defined in (3.10), are Hermitian conjugates of one another. 
Finally, from (3.8) and (3.9), we get 

b+bpRi = eRi (3.13) 

6; 6+ ^2 = eR2 (3.14) 

where 

e= —. 3.15 

We conclude that Ri and i?2 may be considered as eigenfunctions of two SUSY partner 
Hamiltonians [35, 36]. Before using this property to solve equations (3.13) and (3.14) in 
full generality, it is worth reviewing the conditions that should be imposed on Ri and i?2 
to make them physically acceptable. 

As in any bound-state problem for the Dirac equation, the complete relativistic wave- 
function ijj — ^ ^ should be normalizable, which in the present case amounts to the 
condition 

/f ;^(|^i(p)r + |^2(p)f)=l (3.16) 



implying separate convergence of (3.12) for R{p) = R'{p) = Ri{p) and R{p) = R'{p) — 

In deformed quantum mechanics governed by the generahzed commutation relations 
(2.3), it has been shown, however, that normalizability may not be a sufficient condition 
for a wavefunction to be physically acceptable [29]: it should indeed be in the domain of P, 
which physically means that it should have a finite uncertainty in momentum. This leads 
to imposing 



and therefore 



f 

Jo 



dp 2 
p 



dp 2 
p 



f{p) 



Ri{p) + Mp) 



< oo 



Riip) 



< oo 



f 

Jo 



dp 2 

p 



Mp) 



< oo 



(3.17) 
(3.18) 



for the large and small components, respectively. 

In the next section, we will proceed to determine whether such conditions can be fulfilled 
for the minimal e value, i.e., e = or E'^ = 1. 

4 Ground state with E"^ = 1 or e = for given s and 
j values 

From the Hermiticity properties of 6+ and b~ with respect to the scalar product (3.12), it 
results that 



{R,\b;b^\Ri) 

{R2\b;b;\R,) 



{b-Ri\b-Ri) 

{b;R,\b;R,). 



Hence the solutions of equations (3.13) and (3.14) corresponding to e 
solutions of the first-order equations 








(4.1) 
(4.2) 

are necessarily 



(4.3) 
(4.4) 



respectively. Provided they are physically acceptable, these are given by 



(4.5) 
(4.6) 
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where A^f ^ and 7Vi°^ are some normalization coefficients. 

Let us first consider the case of i?f ^ and examine the convergence of the corresponding 
normahzation integral (3.12). Since for p ^ 0, R^^^ behaves as p'', the boundary condition 
Rf\0) = imposes the restriction k > 0. From (3.11), it follows that this can only be 
fulfilled for s = |. Furthermore, for p ^ oo, R^^^ behaves as p~9/f^o^ gQ h^q^i the convergence 
of (3.12) will be ensured provided, in addition, 

i > 4 (-) 

which, for s = |, is equivalent to the condition 

2(3u;j - (3'uj < 2. (4.8) 

Hence is normalizable only for s = \ and small j values satisfying (4.8). 

It is however immediately clear that the convergence of the first integral in (3.18) imposes 
the more stringent condition 

or, explicitly, 

2/3a;(j + 1) + < 2. (4.10) 

We therefore conclude that for R^i^ we must require both s = | and small j values fulfilling 
(4.10). 

Considering next the case of i?2°'* normalizability, from the boundary condition i?2°''(0) = 
we find this time the restriction /c < 0, equivalent to s — —\. Moreover, from the 
behaviour for p — > oo, we get the additional condition q/^q < 1/2, which for s = — | 
amounts to the inequality 2(5ujj — (3'uj < —2. Since for very small values of (3uj and 
such a condition cannot be satisfied by any j value, we conclude that the function ^2°^ is 
never normalizable. 

For those functions Rf^ that are physically acceptable, we must therefore set ^2°^ = 
for the corresponding small component radial momentum wavefunction. In other words, 
we have proved that the set of equations (3.13), (3.14) has well-behaved solutions 
[lif^ 7^ 0,-R2°^ = 0^, corresponding to e = 0, i.e., E = ±1, for s = | and any j value 
satisfying condition (4.10). 

When going back, however, to the coupled radial equations (3.8), (3.9), from which 
the eigenvalue problems (3.13) and (3.14) have been derived, we note that the solutions 
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l^Ri ' 7^ 0, i?2 = Oj are only compatible with E = \ to the exclusion oi E = —\. So, in 
the latter case, although the set of equations (3.13), (3.14) has a well-behaved solution, it 
is unphysical in the sense that it is not a solution of (3.8), (3.9). 

In conclusion, we have shown that as in the case of standard commutation relations [7] , 
physically acceptable wavefunctions of (2.1) with E"^ — 1 can only exist for positive energy 
and 5 = 1, i.e., E — 1 and I — j — \- In the present case of modified commutation 
relations (2.3), there however occurs an additional condition, not present in the standard 
case, namely that the total angular momentum j be restricted to small values satisfying 
equation (4.10). 

In physical terms, this means that there will be neither particles nor antiparticles with 
E = 1 and j values violating condition (4.10), while this remark is only true for antiparticles 
whenever j satisfies (4.10). 

5 Dirac oscillator energy spectrum 

We will now proceed to the solution of equations (3.13), (3.14) (and ultimately to that of 
equations (3.8), (3.9)). In accordance with section 4 results, we have to distinguish between 
the three cases {i) s — ^, small j, (ii) s — ^, large j, and (iii) s — 

5.1 Case | and small j 

It has been shown in section 4 that whenever j satisfies condition (4.10), equation (3.13) 
has a physically acceptable solution Rf\g, k]p) corresponding to the ground state energy 
eo = 0. 

The other eigenvalues corresponding to s = | and the same j value can be found by 
considering ho = bpb~ as the first member of a SUSYQM hierarchy [35, 36] 



j=0 



I = 



0,1,2,... 



(5.1) 



where g^, ki, i = 0, 1, 2, . . . , are some parameters satisfying the conditions 



9o ^ 9 



ki > 




eo = 



(5.2) 
(5.3) 
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and by imposing a SI condition [37, 38] 

bp(gi, ki)bp(gi, h) = ki+i)h~{gi+i, ki+i) + e^+i i = 0, 1, 2, . . . . (5.4) 

Equation (5.4) explicitly reads 

ki+i{ki+i-l) = ki{ki + l) (5.5) 
gi+i{gi+i - Po) = gi{gi + Po) (5.6) 
ei+i = g'i+i(2A;j+i + 1) - 5rj(2A;i - 1) +/3o(A;i+i + /cj). (5.7) 

The solution of the first two equations (5.5) and (5.6), satisfying conditions (5.2) and 
(5.3), is given by 

ki^ k + i g + /3oi. (5.8) 

For such values of Qi and ki, the hierarchy Hamiltonians hi have physically acceptable 
solutions Ri'\gi, ki;p) corresponding to the energies I]}=o^j- 
Using (5.7), we get for the eigenvalues of ho the result 

i 

en = eri{g,k) = ^e^- = An[g + po{k + n)]. (5.9) 
i=o 

Since, in the present case, equation (3.11) yields g — — P{j + |) and k — j + ^, the 
corresponding DO eigenvalues are given by 



E„ — 1 — LU Cn — 4:U!n 



1 - (3u; (^j + + (3ou; (n + j + - 



(5.10) 



In section 6, it will be proved that one can associate to these eigenvalues some physically 
acceptable solutions (^Ri'\ -^2"'') of the coupled equations (3.8), (3.9) for n = 0, 1, 2, . . . or 
n — 1, 2, . . . , according to whether the energy is positive or negative, respectively. 

5.2 Case s = | and large j 

Whenever 

or, in explicit form, 

2Puj{j + l)+p'u;>2 (5.12) 
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there is no physically acceptable wavefunction Ri'\g, k; p) corresponding to e = or = 1. 
We will now proceed to prove that for j values satisfying condition (5.12), there exists a 
ground state of 

ho = b^ig, k)b;ig, k) = - l^f^^ + g{g - (3o)p' + - 2gk - g - (5ok (5.13) 

with e > or > 1. 

For this purpose, it will prove convenient to re-factorize ho. On introducing 

^' = -^ + /3o (5.14) 

equation (5.13) can indeed be rewritten as 

/lo = fe+(^', k)h-{g\ k) + {2k + l){-2g + /?o) (5.15) 

where b^{g', k) are defined as in equation (3.10) but with g replaced by g'. 

If there exists a physically acceptable function Ri'\g', k;p) satisfying the equation 

b;{g\k)Rf^ = (5.16) 
then it will be the ground state wavefunction of ho, corresponding to the eigenvalue 

eo^{2k + l){-2g + Po)- (5.17) 
That this is indeed true follows from section 4 results: R^^ is given by 

where g'Y/^o > 1/2 as a consequence of (5.11) and (5.14), while A; > since s — ^. Hence 
Ri'\g', k;p) is both normahzable and such that its contribution to (p^) converges. 

The remaining eigenvalues of ho can be directly obtained from (5.9) by substituting 
there g' for g, taking (5.14) into account and adding the extra contribution (5.17): 

e„ = H-g + Po{k + l + n)] + {2k + l){-2g + Po)- (5.19) 



On inserting the appropriate values of g and k, the resulting DO eigenvalues read 

1\ „ / 1^ 



-1 + (3uj + -) + PoUJ (n+- 



(5.20) 
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The existence of such eigenvalues is however restricted to that of physically acceptable 
solutions (^Ri^\ -^2""*) of equations (3.8) and (3.9). In section 6, we will show that R'^^ ^ 
contrary to what happens for small j values and £'0 = 1, but that the convergence of its 
contribution to (p^) (i.e., the second condition in (3.18)) imposes a further restriction on 
the allowed g values. This will lead us to distinguish between very large j values for which 
n — 0, 1, 2, . . . are allowed in (5.20), irrespective of the energy sign, and intermediate j 
values for which no bound state does exist. 

5.3 Case s = — | 

Here the negative value of k = —{j + ^) prevents the existence of a normalizable wavefunc- 
tion R^i\g,k]p) corresponding to e = or ii^^ = 1. As in the previous case, however, it 
can be shown by re- factorization that Hq has a ground state with e > or > 1. 
On introducing 

k'^l-k^j + ^ (5.21) 
equation (5.13) can indeed be transformed into 

/lo = h;{g, k%{g, k') + 2(2^ + (3^){j + 1) (5.22) 

where hp{g, k') are defined as in equation (3.10) but with k replaced by k' . 

The ground state wavefunction Rf\g, k';p) of ho, corresponding to the eigenvalue 

eo = 2{2g + (3o){j + l) (5.23) 

is then provided by the solution of the equation 

b;{g,k')Rf^ ^0 (5.24) 

which can be written as 

4°) = N[^)pk'f-(9+Pokm)_ (525) 

Such a function is normalizable since now k' > 0, while 9 = ^ + l3{j + |) satisfies condition 
(4.7). In addition, for (3uj <C 1 and P'lu -C 1, its contribution to (p^) also converges since 
condition (4.9) can now be written as 

2 + 2pu;j>p'cu (5.26) 
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which is fulfilled for any j. 

By proceeding in a similar way as in section 5.2, we find for the remaining eigenvalues 
of ho and the corresponding DO eigenvalues 

en = M9 + /3o(l -k + n)]+ 2{2g + (3,){j + 1) (5.27) 

and 



Ei-l^ u^en = 4a;(n + i + 1) 



1 + (3u + ^) + f^o^ (n + 



(5.28) 



1\ „ / 1^ 

2. 

respectively. In section 6, it will be shown that there exist physically acceptable solutions 
-R2"'') of equations (3.8), (3.9) for any n = 0, 1, 2, . . . , irrespective of the energy sign, 
and that R^'^ ^ 0. 

6 Dirac oscillator radial momentum wavefunctions 

In sections 4 and 5, the ground state wavefunction of Hq — bpb~ has been determined, 
thus providing us with the large component radial momentum wavefunction for the lowest 
state with given s and j values. In this section, we plan to calculate the remaining large 
and small component radial momentum wavefunctions and to examine whether they are 
physically acceptable. 

6.1 Case s = ^ and small j 

The large component radial momentum wavefunctions for the excited states with given s 
and j values correspond to the excited state wavefunctions of ho- According to a well-known 
SUSYQM and SI prescription [35, 36, 37, 38], the latter can be recursively determined from 
the ground state wavefunction (4.5) through the relation 

R^;'\g,k;p)^^^=b;{g,k)R'r-'\9i,h;p) n = l,2,.... (6.1) 

Here wc have specified the {g, A;)-dependence of the operator and wavefunctions and gi — 
g + Po, ki — k + 1, Co = as consequences of (5.8) and (5.9). 

In appendix 2, it is shown that the solution of equation (6.1) is given by 

R^;'\g, k;p) = Ni'-^g, ^)p''+^/-|(«+''+i)pM)(;,) (5.2) 
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where Ni"'\g,k) and P^'^^^z) are a normalization coefficient and a Jacobi polynomial, 
respectively, and we have defined 

^=jX'-\^) "='^-1=^ ^=rf7^ (-!<.<«) (6.3) 
or, conversely, 

1 /2 

, = A(a+i) * = 6 + l (0<P<oc). (6.4) 

For n = 0, wc have shown in section 4 that the corresponding small component radial 
momentum wavefunction R2'\g,k;p) vanishes. For n = 1, 2, . . . , and either En > 1 or 
En < —1, R2''\g,k;p) can be found from Ri^\g,k;p), given in (6.2), by using equation 
(3.9): 

Rt\g,k;p) = -^;^b;{g,k)R'{'\g,k;p). (6.5) 
The result reads (see appendix 2): 

^(")(^, k;p) = Nt\g, fc)p^+^/-^(»+^+i)pf '^)(^) (6.6) 

where 

a = a + l 6 = 6 + 1 (6.7) 

and 

-C/n + -1- 

Note that the n = case may be seen as a special case of (6.6) by defining P^^{^\z) = 0. 

As stated in section 4, (^Rf\g,k;p), R2'\g,k;p) —6^ is a (physically acceptable) so- 
lution of equations (3.8) and (3.9) only for Eq — 1. In contrast, for any n — 1, 2, 
(Rf\g,k]p),Rt\g,k]p)) is a solution of (3.8) and (3.9) for both E^ > I and K < -1. 
It remains to check that such a solution is physically acceptable. 

It is easy to convince oneself that the Jacobi polynomials in (6.2) and (6.6) do not play 
any role in the convergence of (3.16), (3.17), which therefore does not depend on n. To 
analyse this problem it is therefore enough to consider the smallest n value, i.e., n = and 
n = 1, respectively. As a consequence, the convergence is obvious for the large component, 
while for the small one it directly follows from the property i?2^^ ~ -^i^Vp f*^^ P ^ 

Finally, the overall normalization coefficient N["\g,k) can be determined from the 
normalization condition of ip, which amounts to equation (3.16). Use of the Jacobi 
polynomial orthogonality relation [47, 48] then directly leads to N["'\g,k) — [{E^ + 
l)/(2£J„)]^/^A("^(a, 6), where A^'^\a,b) is given by equation (7.8) in section 7. 
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6.2 Remaining cases 

In the two remaining cases corresponding to s — ^ and large j or to s = — |, respectively, 
we can proceed similarly. In recursion relation (6.1), however, we have to replace g by g' 
or k by k' in bp according to which case applies. In contrast, g and k remain unchanged in 
bp when considering equation (6.5). 

As a result, we find that equations (6.2) and (6.6) are still valid provided we appropri- 
ately change the definitions in (6.3), (6.7) and (6.8), as herebelow listed: 

• In the s — ^ and large j case, replace g by g' in (6.2) and define a, b a, b, n and 
Nt\9,k) by 

a = a-l b^b+1 n^n (6.10) 

N!r\s.k) = -aMl±^A,(">(<,',,). (6.11) 

• In the s = — I case, replace k by k' in (6.2) and define a, b d,b, n and N2^\g, k) by 

a — a + 1 b — b — 1 n — n (6.13) 
JVi"'(9.*) = '^f^Nt\!,,k'). (6.14) 

-C/n + l 

As in section 6.1, we still have to check whether the wavefunctions so obtained are 
physically acceptable. Similar arguments as in 6.1 allow us to easily prove this property for 
the large component in both cases and for the small component in the s = —\ case. For 



s 



2 



and large j, however, we note that for p ^ oo, Ri ~ pRf^ ~ p5//3o. As a result, 
although i^g'^'' is normalizable for all g values satisfying condition (5.11), its contribution to 
(p^) is convergent only for 

i<-\ 

or, in explicit form, 

2P(jjj - P'u > 2. (6.16) 

We conclude that physically acceptable solutions of equations (3.8) and (3.9) exist for 
n — 0, 1, 2, ... and either positive or negative energy whenever s = — |ors=| and j 
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assumes very large values satisfying condition (6.16). In the case where s — ^ and j takes 
on some intermediate values in the range 



2-2Puj- P'u < 2Pujj <2 + P'uj 



(6.17) 



no physically acceptable solutions do exist. 

7 Final results 

We may now collect the results obtained in sections 5 and 6 for the energy spectrum and 
wavef unctions. 

Going back to the initial deforming parameters /3, we obtain for the DO complete 
energy spectrum 



Aeon 



1 + (5ujn + (5'uj {n + j + ]^ 



if s = -, j small 



Au{n + i + 1) 
Auj{n + j + 1) 



-1 + (5uj{n + J + 1) + (5'uj 
1 + I3uj{n + j + 1) + p'uj (n + 



if s = 2' ^ ^^^^ ^^^S® 



if s = -- 



(7.1) 



where n runs over all values 0, 1, 2, . . . , except for s = |, j small and negative energy 
where n = 1, 2, . . . only. Small and very large j values are defined by conditions (4.10) 
and (6.16), respectively. No bound states are obtained for intermediate j values satisfying 
condition (6.17). 

Equation (7.1) can be rewritten in a compact form 



4a; 



(7.2) 



by introducing 



e — +1 if s = -, 7 small or if s = — 
2' 2 



= -1 if s = -, j very large. 



(7.3) 
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It can also be recast as 



= 2(^{iV-j + - 



1 



1./ 



l + -i3uj{N-j + -]+ -13' u AT + J + - 



if s = -, j small 



= 2a; (^A^ + j + 

if s = -, J very large 



1 



5 



1 + AT + j + - + -13' AT - j + 



1./ 



1./ 



1 + -/3a; AT + j + - + -/3'a; iV - j + - 



if s = (7.4) 



or 



1 = 2c;[A^ + 2-s-e + (l-2s-e)j]|e+^/?a;[A^ + 2-s-e + (l-2s-e)j] 



+ ^/3'a;[Ar + l + s-(l-2s-e)j]| 



(7.5) 



in terms of the principal quantum number A^ = 2n + Z = 2n + j — s. Here one has to take 
into account that for negative energy, s — \ and a given j value, A^ starts at A^ = j + | 
instead of A?" = j — s as in the remaining cases. 

To the eigenvalues Ensj are associated large and small component wavefunctions defined 
in equations (2.20), (3.1) and (3.5), where the radial momentum wavefunctions can be 
written as 



1/2 



— —ea 



Ensj 1 



1/2 



A(")(a,6)p^+5/-iM+i)pf'' 



(7.6) 
(7.7) 



Here 



_ [ 2l3'o-^\a + b + 2n + l)n\T{a + b + n+l) \ 
^ " V r(a + n + l)r(6 + n + l) ) 



1/2 



a — 



(3qu 
a — a + e 



1-I3u 

6 = 6 + 2s 



2 + «(2J + 1) 



n — n — s 



nsj 



\E. 



(7.8) 

(7.9) 
(7.10) 
(7.11) 



nsj I 



and n runs over 0, 1, 2, ... , except for s = |, e = +1 and a — —1, where n = 1, 2, 
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For j3uj, f3'uj — > 0, it is obvious that all j values fall into the "small" category, as defined 
by equation (4.10), so that equations (7.4) and (7.5) become 



2.(iV-i + i) i£. = i 

2c^(jV + i + |) ifs = -i (7.12) 



and 



El^^ - 1 = 2uj[N + 1 - s{2j + 1)]. (7.13) 

Hence the deformed DO energy spectrum smoothly goes over to the nondeformed one [5, 7]. 
In the latter, there is an infinite (resp. a finite) number of degenerate states with s — \ 
(resp. s = — |), corresponding to a given N — j (resp. N + j') value [5]. As a consequence, 
the conventional DO has an so(4) ® so(3, 1) dynamical symmetry Lie algebra [8]. 

For nonvanishing (5uj and (5'uj values, this degeneracy scheme is completely spoilt so that 
the conventional DO symmetry Lie algebra gets broken. 

8 Conclusion 

In the present paper, we have shown that the DO problem remains exactly solvable in 
the momentum representation after modifying the canonical commutation relations so that 
there appear isotropic nonzero minimal uncertainties in the position coordinates. The 
factorization method (or, equivalently, SUSYQM and SI techniques) has proved very con- 
venient for deriving both the energy spectrum and the corresponding wavefunctions of this 
generalized DO problem. 

The outcome of our study shows some resemblances to the conventional results. 

One of them is that there is no symmetry under exchange of s = | with s — —\. Such 
a dissymetry is related to the specific substitution P ^ P — icoXp that has been used 
in constructing the DO Hamiltonian (2.1) and which is only one of the two possibilities 
P — > P =F iujXP at our disposal. On considering indeed equations (3.8) and (3.9) for 
uj replaced by — cu, it results from the properties bp{—u!,s) = —bj{uj,—s) that they get 
equivalent to the same where s is changed into —s, and E, Ri, R2 are replaced by —E, 
—R2, Ri, respectively. It should be stressed that in the conventional case, the u — > —uj 
transformation has been considered in connection with supersymmetry [11] or, equivalently, 
with the particle antiparticle transformation [14]. 
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Another similarity, related to supcrsymmctry too, lies in the absence of negative-energy 
states with vanishing energy cqsj = (i.e., Eosj = — 1). 

There are also some important differences between the deformed and conventional DO 
problems. 

When modifying the canonical commutation relations, there appear, as expected, some 
additional terms in the explicit expression for E^^j — l, which becomes quadratic in n instead 
of linear. As a consequence, the degeneracy pattern observed for the conventional DO [5, 7] 
gets completely spoilt, so that the corresponding dynamical symmetry [8] is broken. 

More unexpectedly, deformation leads to a new and interesting feature in the energy 
spectrum: for s = i, there appears a difference in behaviour between distinct j values. 
For small ones defined by condition (4.10), the ground state has vanishing energy e^sj = 
{Eqsj = 1) and SUSY is unbroken as in the conventional case. For very large j values such 
that condition (6.16) is satisfied, the ground state acquires a nonvanishing energy e^sj > 
{Eq^j > 1), as in the s = — | case, so that SUSY is broken. Finally, for intermediate j 
values in the range (6.17), the DO has no bound state . For deriving all these properties, it 
has been essential to impose not only that the (relativistic) wavefunctions are normalizable, 
but also that they lead to a finite uncertainty in P, a supplementary condition specific of 
the type of deformation considered here. 

There remains an interesting open question in connection with the use of the deformed 
commutation relations (2.3) in the DO relativistic problem. Since the conventional DO is 
Lorentz covariant [6], one may wonder whether such a property remains true after defor- 
mation. Although it has been known for a long time that quantized space-time may be 
compatible with Lorentz invariance [49], it is clear that the formalism based on relations 
(2.3) is not. As already noted before [50], a boost may indeed squeeze any "minimal" 
length, such as A.Xq, as much as one likes. From a mathematical viewpoint, the question is 
therefore whether the relations (2.3) can be associated with some deformed Lorentz algebra 
(as occurring, e.g., in the Minkowski non-commutative space-time [51]) and, in such a 
case, what would be the transformation properties of the deformed DO. From a physical 
viewpoint, this is related to a question of considerable interest, which has been most debated 
in the recent quantum-gravity literature, namely the possibility of Planck-scale departures 
from Lorentz symmetry and its experimental testing (for a recent review see, e.g., [52]). 
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Appendix 1. The operators in momentum repre- 
sentation 



The purpose of this appendix is to prove equations (2.16) and (2.17), as well as equation 
(2.27). 

On inserting (2.8) in definition (2.13) of B^, the latter become 



B = a ■ p^uj 



d 



d 



1 + l3p^)(T ■ — + I3\a ■ p) p ■ ^ + • P 



dp \ dp 

To proceed further, it is convenient to apply the well-known relation 

{a Oi){a O2) = Oi O2 + i<T (Oi x O2) 



(Al.l) 



(A1.2) 



valid for any two vector operators Oi, O2, whose components commute with those of a. 
It leads to the following results: 

d 



[(T -p) \a 



d 



dp. 



a.- (a-p) 



or, alternatively. 



d . d\ 
dp \ dp J dp 

d f d \ d 

— - p + ia-l — xp\= p— + CT- L + 3 
dp \dp J dp 



d a L 



(A1.3) 
(A1.4) 



dp p 
d 1 (T- L + 3 d (T-L + 2 

dpp p dp p 



(A1.5) 
(A1.6) 



where use has been made of (7p, defined in (2.18). Property (2.19) then allows one to write 
the operator cr ■ d/dp, appearing in equation (Al.l), as 



d 



d a L 



d ^ (T-L + 2 ' 



(A1.7) 



dp ^ \dp p ) \dp p 

Furthermore, the operator (cr ■ p){p ■ d/dp), also appearing in equation (Al.l), can be 
expressed as 



(fT -p) {p~] ^app^-^^p'^-^ap (A1.8) 



d_' 
dp, 



dp dp ^ 



as a consequence of the relation 



d_ _ d_ 
dp ^ ^ dp 



(A1.9) 
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which is easily derived from equation (2.18). 

On combining (Al.l) with (2.18), (Al.7) and (Al.8), it is then straightforward to get 
equations (2.16) and (2.17). 

Let us next consider equation (2.27). From (Al.3) and (Al.4), the operator cr • L can 
be expressed as 



d\ d 
cr ■ — \ pap - p- 3 



dp 



dp 



from which it follows that 



These relations lead to the equation 

{(7p, a- ■ L} 

equivalent to (2.27), which is therefore proved 



d 

op 



— ScTp — — 2(7p 



(ALIO) 



d f d \ f d\ d 

ap{a ■ L) = pap— - p I a ■ — 1 {(T ■ L)ap = \ cr ■ — \ p - pap— - ?,ap. (ALU) 



(A1.12) 



Appendix 2. Calculation of radial momentum wave- 
functions 



The purpose of this appendix is to give some details on the determination of the radial 
momentum wavefunctions carried out in section 6 . 

Let us first consider equation (6.2) for s = | and small j. On using (6.3), it is clear that 
for n = it gives back the ground state wavefunction (4.5), as it should be. Furthermore, 
insertion of (6.2) (where n ^ n — 1, g ^ gi and k ki) in the right-hand side of (6.1) 
leads to the relation 



Rt\g,k-p) 
Nt'\guh) 



'^(gl,fcl) b+i r-Ua+b+l) 



1- z^)^ + a-h+{a + h + 2)z 



dz 



2nAr(" '\gi,ki) ,^1 .-lia+b+i) p{a,, 



(A2.1) 
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where, in the last step, we have used equation (1.8.7) of Ref. [47]. This completes the proof 
that equations (6.2) and (6.3) provide the sohition of (6.1). 

Let us now turn ourselves to equation (6.6). On combining (6.5) with (6.2), we obtain 



En + 1 



(JjZ 



Er, + 1 



^ 2/3oa;(a + 6 + n + l)N^{g, A:) ^fe+3^_ 1(^+^+3) ^^a+i,j,+i)^^^ ^^2.2) 

demonstrating the correctness of equations (6.6) - (6.8). Here, in the last step, we have 
taken advantage of equation (10.8.17) of Ref. [48]. 

In the s = I and large j case, equation (6.5) for the small component radial momentum 
wavefunction becomes 



^ uN^{g',k) 
En + 1 



,d , / 1\ b + l 

d 



^ 2(3ou;N-{g',k) .-i(a+.+i) 



2(5^uj{a + n)N^{g',k) 



(A2.3) 



En+l 

on combining various properties of Jacobi polynomials [48]. Equation (A2.3) agrees with 
equations (6.6), (6.10) and (6.11). 
Finally, in the s — —\ case, we get 



^ ujN^{g,k') 

En + l 

2ujN^{g,k'] 



\d ^ ( 1\ h-\ 



En + \ 

^ 2^(& + n)7Vf(g,A;0 ^,_i^_i(„^,^,)^^„^,^,_,)^^^ 

En + 1 



/-5(«+''+l)p^«'^)(^) 
p(a,6)(^) 



(A2.4) 



on combining various properties of Jacobi polynomials again [48]. Hence equations (6.6), 
(6.13) and (6.14) are valid. 
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